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ABSTRACT 

We compare non-locality of interactions between different scales in hydrodynamic (HD) turbu- 
lence and niagnetohydrodynamic (MHD) turbulence in a strongly magnetized medium. We use 3- 
dimensional incompressible direct numerical simulations to evaluate non-locality of interactions. Our 
results show that non-locality in MHD turbulence is much more pronounced than that in HD turbu- 
lence. Roughly speaking, non-local interactions count for more than 10% of total interactions in our 
MHD simulation on a grid of 512"^ points. However, there is no evidence that non-local interactions 



are important in our HD simulation with the same numerical resolution, 
non-locality affects energy spectrum. 

Subject headings: ISM:general — MHD — turbulence 
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1. INTRODUCTION 

Turbulence is commonly observed in astrophysical flu- 
ids and in many cases such turbulence is accompanied by 
a strong magnetic field, which has a large impact on the 
dynamics of the turbulent cascade. Since turbulence in- 
fluences many astrophysical processes (e.g. transport of 
mass and angular momentum, star formation, fragmen- 
tation of molecular clouds, heat and cosmic ray trans- 
port, magnetic reconnection, etc.), understanding scaling 
properties of niagnetohydrodynamic (MHD) turbulence 
is essential for theoretical astrophysics. For this reason, 
rich literature exists regarding scaling relations of MHD 
turbulence (See Goldreich & Sridhar 1995; Biskamp 
2003 and references therein; see also Cho & Vishniac 
2000b; Maron & Goldreich 2001; Miiller, Biskamp, & 
Grappin 2003; Miiller & Grappin 2005; Boldyrev 2005, 
2006; Beresnyak & Lazarian 2006; Mason, Cattaneo 
& Boldyrev 2006; Gogoberidze 2007; Matthaeus et al. 
2008). 

In hydrodynamic (HD) turbulence, energy cascades 
down to smaller scales. Kinetic energy contained in 
an "eddy" is transferred to smaller eddies by shearing 
motions of other eddies (see, for example, Frisch 1995). 
Most theories on turbulence assume locality of interac- 
tions, which means interactions between similar size ed- 
dies dominate in such energy cascade. In Fourier space, 
this means that a Fourier mode at a wavenumber fc = |k| , 
where k is the wavevector, interacts mainly with other 
modes having similar wavenumbers and transfers its en- 
ergy to modes that have larger wavenumbers. Recently 
many researchers have investigated locality in HD tur- 
bulence (Mininni, Alexakis, & Pouquet 2008; Alexakis, 
Mininni, & Pouquet 2007; see also Verma et al. 2005). 

In MHD turbulence with a strong mean field (Bq), lo- 
cality is also generally assumed. However, in the MHD 
case, the nature of energy cascade is slightly different. 
In the incompressible limit, any magnetic perturbation 
propagates along the magnetic field line. To the first or- 
der, the speed of propagation is constant and equal to 
the Alfven speed Va = -Bo/^/47rp, where p is the den- 
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sity. Since wave packets are moving along the magnetic 
field line, there are two possible directions for propa- 
gation. If all the wave packets are moving in one di- 
rection, then they are stable to nonlinear order (Parker 
1979). Therefore, in order to initiate turbulence, there 
must be opposite-traveling wave packets and the energy 
cascade occurs only when they collide. Therefore, in the 
MHD case, locality means that a wave packet (or "eddy") 
transfers energy to smaller scale wave packets by shear- 
ing motions of opposite-traveling wave packets of similar 
size. 

There have been some discussions about non-locality in 
MHD turbulence with a strong mean fieldr] For example, 
Alexakis (2007) theoretically studied nonlocal model of 
MHD turbulence. In their inspiring work, Beresnyak & 
Lazarian (2010) numerically studied non-locality and ar- 
gued that "MHD turbulence is fairly non-local, at least 
less local than hydrodynamic turbulence" (see also Beres- 
nyak & Lazarian 2009). They claimed that " a) the lack 
of visible bottleneck effect in MHD turbulence, while it 
is clearly present in hydro turbulence, and b) the depen- 
dence of kinetic and magnetic spectra on driving" sup- 
port this idea. Teaca et al. (2009) calculated anisotropic 
energy transfer in Fourier space. But they did not pay 
much attention to the locality issue. 

In this paper, we quantitatively evaluate non-locality 
of HD and MHD turbulence and present a direct evidence 
that non-locality is clearly present in MHD turbulence. 
We consider only balanced strong MHD turbulence. Here 
balanced MHD turbulence means that amplitudes of two 
opposite-traveling wave packets are almost equal. In §2, 
we describe our numerical setup. In §3, we present our 
results. In §4, we briefly discuss how non-locality affects 
energy spectrum and give summary. 



2. SIMULATIONS 
We solve the incompressible HD equation, 

dt-v ^ -(V X v) X V + i^V^v + f - VP', 



(1) 



^ When the mean field is weak or zero, turbulence structure 
is very different (see for example Cho et al. 2009). There are 
many discussions about non-locality in this regime (see for example 
Alexakis, Mininni, & Pouquet 2005b; Lessinnes, Verma, & Carati 
2008; Yousef et al. 2009; Aluie & Eyink 2010). 
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Fig. 1. — (a) Time evolution of v'^ and h^ . (b) Compensated 
spectra at t = 12. 

and the incompressible MHD equations, 
dt'v = -(Vxv)xv + (VxB)xB + i^V^v+f-VP', (2) 



9tB = B • Vv - V VB + T^V^B, 



(3) 



in a periodic box of size 27r, where f is a random forcing 
term with unit correlation time, P' = P + u^/2, P is 
pressure, v is the velocity, and B is the magnetic field 
divided by (47r/9)^/^. Thus the field B is, in fact, the 
Alfvenic velocity. The velocity and the magnetic fields 
are divergence-free: V-v ~ V-B = 0. The peak of energy 
injection occurs at fc^ w 2.5, so the energy injection scale 
is L ~ 2.5. The amplitudes of the forcing components 
are tuned to ensure v ~ 1. 

In the MHD simulation, the magnetic field consists of 
the uniform background field and a fluctuating field: B = 
Bo + b. The Alfven velocity of the uniform background 
field, Bq, is set to 0.8. At i = 0, the magnetic field 
has only the uniform component. We consider only the 
case where viscosity is equal to magnetic diffusivity: v = 
■q. Details of the code can be found in Cho & Vishniac 
(2000ab). 

Figure [ija) shows time evolution of kinetic and mag- 
netic energy densities. Figure flVb) shows energy spectra 
at t=12. The kinetic spectrunifor HD Run (solid curve) 
is consistent with the Kolmogorov spectrum {E{k) ex 
^-5/3-j £qj. y. g (2,15). But it shows a moderate in- 
crease of the slope for k > 15. The kinetic spectrum 
for MHD Run (dashed curve) is also consistent with the 
Kolmogorov one. However, the magnetic spectrum (dot- 
ted curve) is slightly shallower than the Kolmogorov one. 
Therefore, the spectrum of b"^ +v^ (not shown) is slightly 
shallower than the Kolmogorov one. 

3. RESULTS 



3.1. Shell-to-Shell Energy Transfer 

We can rewrite the MHD equations in Eqs. |2|) and 
([3| using the Elsasser variables, Z+ = v + B antiZ^ = 



B: 



StZ+ = -Z" ■ VZ+ 

a,z^ = -z+ -vz" 






f - VP, 
- f - VP 



(4) 
(5) 



The Elsasser variables denote amplitudes of two 
opposite-traveling waves along the magnetic field line. 
The nonlinear term of Z+, for example, states that en- 
ergy transfer between Z+ modes is mediated by Z~ 
modes. 

In Fourier space, the non- linear term in Eq. Q, for 
example, becomes 



N+ = -ik • 



E 

p-fq=k 



ZpZj 



and the time derivative of (l/2)|Zj^p is given by 



K-dt'LV 



= Z+.N+^ 



/fc- 



\K\\ 



(6) 



(7) 



where '*' denotes complex conjugate and we dropped the 
forcing term because its role is limited in the inertial 
range. Energy transfer occurs only between Zj^ and Z+, 

via shearing motions provided by Z" modes. Without 
Tr modes, Z+ modes alone do not interact each other. 

If interactions are local in Fourier space, we will have 
p ^ q r^ k. Since it is difficult to check the locality using 
individual triad interactions in Fourier space, we investi- 
gate shell-to-shell interactions. That is, we consider col- 
lective interactions in Fourier space between Z+ modes 
in a unit shell of radius k (hereinafter, "fc-shell") and Z+ 
modes in a unit shell of radius q ( "(/-shell" ) by the help 
of Z^ modes in a unit shell of radius p ("p-sheU")rl 

We first consider shell-to-shell interactions of Z+ 
modes mediated by all Z~ modes. Contour diagrams in 
Figure l2ja) and (d) show the shell-to-shell energy trans- 
fer rate: 



T{k,p = all, q) = ^ T{k,p, q) = 



(8) 
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where p = k' — q' and Pmax is the largest wavenumber, 
for MHD (Figure 2jd)) and a similar expression for HD 
(Figure [2][a))Q The energy transfer rate T{k,p = all, q) 
here is similar to T2{K, Q) in Alexakis, Mininni, & Pou- 
quet (2005a) or Mininni et al. (2008). The contour di- 
agrams are exactly anti-symmetric with respect to the 
k = q line. The overall shape of the contour diagram 
for HD is consistent with earlier findings (Alexakis et al. 
2005a; Mininni et al. 2008). The value of T(fc,p = all, q) 
is positive on the upper- left half, which means that, when 
q < k, the Z"*" modes in fc-shell gains energy from the Z^ 
modes in (7-shell by the help of all Z~ modes that satisfy 
p + q = k. This result is consistent with the concept of 

^ In this paper, wavenumbers p, q, and k refer to those described 
here (and those appear in Eqs. Kh, (Tfjl and l|8|l). 
" In the HD case, Z+ = Z" = v. 
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Fig. 2. — (a) & (d): Energy transfer from g-shell to fc-shell by mediation of all p modes that satisfy p = k — q. (b) & (e): Energy transfer 
from g-shell to fc-shell by mediation of p modes that satisfy 1.5 < |p| < 3.5. This is equal to the Energy transfer from q-shell to fc-shell 
mediated by the outer scale motions, (c) & (f): T{ko,*,q) vs. q for fco = 6, 16, and 32. The values for fco = 16 and 32 are offset by 0.015 
and 0.03, respectively. The solid lines are for T(fco,p = all, g) and the dotted lines for T{ko,p = 2&3,g). In the MHD case, the solid and 
the dotted lines look very similar, which might mean that the outer scale counts for substantial fraction of the energy transfer rate. In 
contour diagrams the linear color bars range from -0.005 to 0.005. Snapshot at t = 12. 



energy cascade: energy cascades down to smaller scales. 
Note that the values of T{k,p — all,q) are very close 
to zero except near the k = q line. Does it mean that 
locality is a good approximation? 

Note that, when the outer scale of turbulence provides 
strong shearing motions, T{k,p = aU,q) has positive 
peaks at q = k — ki, and negative peaks at q = k + k^. 
In our case, k^ ^ 2.5. Therefore, shearing motions of 
the outer scale can also produce diagrams similar to Fig- 
ure [2la) and (d). Indeed, when we plot 



r(fc,p = 2&3,9) = ^r(A;,p,q), 



(9) 



which is similar to T(k,p — all, q) except the fact that we 
do the summation from p — 2 to p ~ 3, the contour dia- 
grams show similar features (Figure [2jb) and (c)). This 
result is consistent with earlier results for HD turbulence 
(e.g. Alexakis et al. 2005a). Note that, in Figurel2ja) and 
(d), the width of the contour lines near the k = q line 
is narrower in the MHD case than in the HD case. This 
might mean that the effect of the outer scale is stronger 



in the MHD case than in the HD case. 

Figure [2lc) and (f ) (right panels) show the values of 
T{ko,p = all,q) (solid) and T(fco,P = 2&3, g) (dotted) 
for select values of fep. We take fcg — 6, 16, and 32. The 
values for fcg = 16 and 32 are offset by 0.015 and 0.03, 
respectively, for clarity. In the HD case (upper panel), 
T(fco,p — all,q) and T{ko,p = 2&:3,(7) look different. 
However, in the MHD case (lower panel), T{ko,p = all, q) 
and T{ko,p = 2&3, g) look very similar, which might 
mean that the outer scale does play important roles in 
MHD energy cascade. 

In order to evaluate the role of the outer scale in shcU- 
to-shell energy transfer, we calculate the ratio 



fe-i 



fc-i 



Y, T{k,p = 2&3, q)/ J2 T{k,p = all, q). 



(10) 



9=0 



q=a 



where qmax is the largest wavenumber. Figurepla) shows 
the ratios for HD and MHD. In the HD case (solid curve), 
the ratio is less than 0.5 for most values of k, which is 
consistent with Mininni et al. (2008). However, in the 
MHD case (dotted curve) the values are > 0.5 for most 
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Fig. 3. — Ratio of local to non-local interactions, (a) The ratio of energy transfer from all g-shells with g < A; to fc-shell mediated by 
the outer scale motions (p = 2 and 3) to that mediated by all scale motions (all p-shells). (b) The ratio of energy transfer from all g-shells 
with q < kjnin = A;/\/2 to all fc-shells between fc/\/2 and \/2fc (we call "fc-band"), mediated by the outer scale motions (p = 2 and 3) to 
that mediated by all scale motions (all p-shells). Snapshot at t = 12. 

values of k, which means that non-local interactions are 
indeed important for shell-to-shell energy transfer in the 
MHD case. 

However, it is very important to note that the result 
in Figure p[a) does not mean that non-local interactions 
are as strong as local interactions in MHD cascade. The 
result in Figure IsFa) is only for a single shell. 

In order to evaluate non-locality, we would better con- 
sider the effect of p-shells on a band of fc-shells between 
kmin = k/a and kmax — ctk, where a is a constant. In 
this paper, we take a = v 2. The motivation for consid- 
ering this quantity is that Fourier modes in {k/\/2, \/2k) 
can define "eddies" on a scale I ^ 1/fc. In FigurelsFb) we 
plot the ratio similar to that in Eq. (10), but expressed 
in terms of 



that non-locality is indeed present in MHD turbulencejj 

3.2. More on Non-locality of MHD Turbulence 

Figure Islb) shows that the p — 2 and p = 3 shells 
contribute more than 10% of the total energy flux. Then, 
which p-shell provides the strongest contribution to a fc- 
band? In other words, what is the most shear-providing 
shell for a band of fc-shells between k/^/2 and ^/2k7 To 
see this, we calculate the following quantity: 



24 



24 



^T{ki,and,P,q) ^^ ^ T{k',p,q), 



(12) 
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(11) 



k' = k„ 



where k,nin == fc/v2 and kmax 



V2fc. The summation 
for q is done from to kmin — 1- The ratio for MHD 
(dotted) is non-negligible and still substantially larger 
than that for HD (solid). Therefore, we can conclude 



which is equal to the total energy transferred from all 
g-shells with g < 24 to the fc-band between fc = 25 and 
fc = 50 by the shearing action of a p-shell. Figure |4] shows 
that each p-shell provides a similar contribution in the 
HD case (solid line). Therefore, non- locality does not 

® We note that the ratio for MHD gradually decreases as k 
increases. Although it is not very clear at this moment whether 
it will continue to drop when we have a very long inertial range, 
it is likely that the ratio will continue to drop and the non-local 
effects of the outer scale will ultimately vanish on very small scales. 
Nevertheless, non-locality is an important characteristic of MHD 
turbulence near the outer scale. 
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Fig. 4. — The amount of energy transferred to a fc-band (25 < 
k < 50) from all tj-shells having wavenumbers smaller than 25 (i.e. 
q < 24) as a function of p, which is the radius of a p-shell in Fourier 
space. Solid line is for HD turbulence and dotted line for the MHD 
case. In HD, non-locality is not conspicuous. In MHD, we can see 
a sharp peak at p = 2. Snapshot at t = 12. 



seem to be important in HD turbulence. However, the 
p ^ 2 shell contributes most in the MHD case (dotted 
line) . This is another piece of evidence that non-locality 
is clearly present in MHD turbulence. 

3.3. Non-local energy transfer 

So far, we have discussed non-local influence of the 
outer scale eddies. Now, it is time to clarify the meaning 
of non-locality. The nonlinear term for 9tZ+, for exam- 
ple, is — Z^ • VZ+. Since the nonlinear term contains 
both Z and Z+, non-locality has two meanings: 

1. Non-local effects of Z~, and 

2. Non-local effects of Z+. 

Since Z~ modes do not lose or gain energy, the former 
type of non-locality does not involve energy transfer be- 
tween the outer scale and small scales. In fact, when 
this kind of non-locality is present, energy transfer be- 
tween adjacent shells is enhanced. Non-locality we have 
discussed so far is this type of non-locality. 

When the latter type of non- locality is present, there is 
direct energy transfer between different scales. In order 
to evaluate the energy transfer rate from the outer scale 
to a fc-band between fc^m = fc/"\/2 and k^ax = V^k, we 
calculate the ratio 



T{kband,P = all, q = 2&3) 



Z^n'=0 ' [kband,P 



all, q') 



(13) 



Figure IsFa) shows the ratios for HD and MHD. We can 
see that the ratios for this type of non- locality are smaller 
than those for the former type of non-locality (see Fig- 
ure Islb)). The ratio for the MHD case is higher than 
that for the HD case. 



50 



T{kband,P^all,q)^ 2_^ Tik',p = all,q) (14) 

fc'=25 

for HD and MHD. The values of T{kband,P — all, 9 = 2) 
and T{ki,and,P — all,q — 3) are not particularly larger 
than other values. It is clear from the figure that energy 
transfer from the outer scale (g = 2 and 3) to the fc-band 
(25 < fc < 50) is small. Therefore, non-local energy 
transfer from the outer scale to small scales may not 
be an important characteristic for both the HD and the 
MHD cases. 

4. DISCUSSION AND SUMMARY 

When shearing motions of the outer scale eddies influ- 
ence energy transfer of inertial range eddies, energy spec- 
trum becomes flatter than the Kolmogorov one. Suppose 
that the shearing motions of the outer scale eddies com- 
pletely dominate energy cascade. In this case, from 



Zi /teas ^ Zi /{L/vl) (x Zi = constant, 



(15) 



we can easily show that energy spectrum is E{k) ex fc~^, 
where Zi is an Elsasser variable at scale I, L the outer 
scale, and vl the rms velocity at the outer scale (see 
Equation (1) of Cho, Lazarian, & Vishniac 2003). If 
the shearing motions of the outer scale eddies do not 
completely dominate, we will have a spectrum between 
fc-i and fc-5/3. 

Indeed, in Figure fl] we observe that Ey{k) + Eb{k) 
(hence spectrum of Z^ or Z~) in MHD is flatter than 
the Kolmogorov spectrum. This is consistent with earlier 
numerical results (see, for example, Maron & Goldreich 
2001; Miiller et al. 2003). 

In summary, we have found the following results. 

1. We have developed a quantitative method to mea- 
sure non- locality (see Figure IsFb) or Figure El). 



2. Our numerical calculations show non-locality is 
more pronounced in MHD turbulence than in HD 
turbulence. This result confirms an earlier finding 
by Beresnyak & Lazarian (2010). 

3. There are two forms of non-locality in MHD tur- 
bulence: non-local influence of shearing motions, 
which does not involve energy transfer between dif- 
ferent scales, and non -local energy transfer between 
different scales (j |3.3[ ). In MHD, the former type of 
non-locality (i.e. non-local influence of outer scale 
shearing motions) is more important. It is not clear 
whether the latter type is important. 

4. In MHD, non-locality is not negligible so that it 
might affect dynamics of turbulent cascade. 
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Fig. 5. — Non-local energy transfer, (a) Ratios of local to non-local direct energy transfer rates. The quantity shown here, T{k,p = 
all,q = 2&3)/T(fc, aii, < g < kmin), is the ratio of non-local energy transfer rate to total energy transfer rate. The numerator is the 
rate from q = 2 and q = 3 shells to fc-bands between kmin = k/%/2 and kmax = \f2k and the denominator is the total energy transfer 
rate from (/-shells between and kmin — 1 to the same fc-bands. The MHD case shows stronger non-locality, (b) The amount of energy 
transferred from a g-shell to a fc-band (25 < fc < 50) by the mediation of all p-shells. Solid line is for HD turbulence and dotted line for 
MHD turbulence. In both cases, non-locality is not conspicuous. Snapshot at t = 12. 
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